Effective Actions for the SU(2) Confinement— Deconfinement Phase Transition 



in 
o 
o 

(N 



> 
m 

o 

(N 
O 

o 

in 

<D 
43 



X 



Thomas Heinzl* 
School of Mathematics and Statistics 
University of Plymouth 
Drake Circus, Plymouth, PL4 8AA, United Kingdom 

Tobias Kaestner'f and Andreas Wipf* 
Theoretisch-Physikalisches Institut, 
Friedrich-Schiller-Universitat Jena 
Max-Wien-Platz 1, 07743 Jena, Germany 
(Dated: February 1, 2008) 

We compare different Polyakov loop actions yielding effective descriptions of finite-temperature 
SU(2) Yang-Mills theory on the lattice. The actions are motivated by a simultaneous strong- 
coupling and character expansion obeying center symmetry and include both Ising and Ginzburg- 
Landau type models. To keep things simple we limit ourselves to nearest-neighbor interactions. 
Some truncations involving the most relevant characters are studied within a novel mean-field ap- 
proximation. Using inverse Monte-Carlo techniques based on exact geometrical Schwinger-Dyson 
equations we determine the effective couplings of the Polyakov loop actions. Monte-Carlo simula- 
tions of these actions reveal that the mean-field analysis is a fairly good guide to the physics involved. 
Our Polyakov loop actions reproduce standard Yang-Mills observables well up to limitations due to 
the nearest-neighbor approximation. 
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I. INTRODUCTION 

The finite-temperature confincment-deconfinement 
phase transition in SU(2) Yang-Mills theory originally 
conjectured by Polyakov [1] and Susskind [2] is by now 
fairly well established. The order parameter is the 
Polyakov loop, 



1 / f^ T \ 

L x = -trT exp\i J dTA°{x,r)j , 



(1) 



a traced Wilson line that winds around the periodic Eu- 
clidean time direction parameterized by t, < r < Pt 
where (3t — 1/T is the inverse temperature. In the con- 
fined phase the expectation value (L) is zero, while it 
becomes nonvanishing in the broken, deconhned phase. 
The Polyakov loop transforms nontrivially under the cen- 
ter symmetry, 



L x ^zL x , z = ±1 G Z(2) 



(2) 



Thus, above the critical temperature, T = T c , this sym- 
metry becomes spontaneously broken. Lattice calcula- 
tions have shown beyond any doubt that the phase tran- 
sition is second order with the critical exponents being 
those of the 3d Ising model [3-6]. This is in accordance 
with the Svetitsky-Yaffe conjecture [7, 8] which states 
in particular that SU(2) Yang-Mills theory (in Ad) is 



in the universality class of a Z(2) spin model (in 3d) 
with short-range interactions. Hence, it should be possi- 
ble to describe the confinement-deconfinement transition 
by an effective theory formulated solely in terms of the 
Polyakov loop L x . The most general ansatz is given by a 
center-symmetric effective Polyakov loop action (PLA) 
of the form [9] 

s pl [l] = j2nLl]+J2 L * K ^ L y 

x xy 

+ ^ L x L y K^j uv L u L v + . . . (3) 

xyuv 

There is a potential term V , a power series in L x living 
on single lattice sites x, plus hopping terms with kernels 
j( (2n) connecting more and more lattice sites x,y, . . .. It 
is important to note that L is a dimensionless and com- 
pact variable, L x G [—1,1]. Thus, in principle, one is 
confronted with a proliferation of possible operators that 
may appear in the PLA (3). One simplification arises 
due to the Svetitsky-Yaffe conjecture implying that the 
kernels K ly2n%> should be short-ranged. Hence, upon ex- 
panding like for instance, 



K, 



(2) _ 



xy 



y,x-\-r 



(4) 
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one expects that the first few terms with small r = \x — y\ 
will dominate, i.e. will have the largest couplings A r . To 
check this expectation one needs a reliable method to 
calculate the kernels X( 2 ") or, equivalently, the coupling 
parameters inherent in them. A particularly suited ap- 
proach is introduced in the following section. 
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II. INVERSE MONTE-CARLO METHOD 



given by the product of temporal links, 



Inverse Monte-Carlo (IMC) is a numerical method to 
determine effective actions [10, 11]. The latter are gener- 
ically defined via 



y x = I] tf-.t.o e su(2) . 



(7) 



t=l 



cxp(-S* cff [X]) = J&U8(X-X[U])exp(-S[U\) , 



(5) 



This SU(2) matrix may be diagonalized whereupon it 
can be written as 



where the U's represent some 'microscopic' degrees of 
freedom and the X's the effective 'macroscopic' ones. In 
the spirit of Wilson's renormalization group, these are 
obtained by integrating out the U's in favor of the X's. 
It is important to distinguish this 'Wilsonian' notion of 
an effective action from the 1 -particle-irreducible (1PI) 
effective action which will later on be employed as well 
(cf. the recent remarks in [12]). 

Of course, the problem with (5) is to do the integra- 
tion nonperturbatively which in general is not possible. 
In this case, one has to resort to choosing an ansatz like 
(3) as dictated by symmetry and dimensional counting 
or to do the integration numerically. The huge num- 
ber of degrees of freedom involved clearly suggests to use 
Monte-Carlo (MC) methods. However, this amounts to 
calculating expectation values rather than integrals like 
(5). Hence, one needs a recipe to get from expectation 
values to effective actions. This is exactly what IMC is 
supposed to do (see Fig. 1). 

integrate out 

action S[U] * off. action SWrfA] 

U 



m _ ( exp(i9 x ) 






exp(-z6> x ) 



with — 7r < 8 X < ir. This representation immediately 
yields the trace (divided by two), 



itr<p a 



cos Q a 



(8) 



which contains all the gauge invariant information con- 
tained in the group variable ty x . As an aside, we remark 
that this peculiar feature will no longer be true for higher 
SU(N) groups [13, 14]. In this more general case, traces 
of N — 1 different powers of ^ are required. 

With (5) the action (6) leads to the partition function 



Z = J ®L exp(-S PL [L]) 



(9) 



where the integration is performed with the reduced Haar 
measure h of SU(2), 



@L = Y[dh{L x ), dh{u) = -^/l-u 2 du . (10) 



MC 



conf's {U} 



IMC 



calculate 



X = X[U] 

FIG. 1: Illustration of the IMC procedure. 

Our particular IMC method is based on the Schwinger- 
Dyson equations that must hold in the effective theory 
once a particular ansatz is chosen. In our case, the macro- 
scopic degrees are given by the Polyakov loops L x dis- 
tributed according to the PLA (3) which we rewrite as 



Since SpL depends on the Polyakov loop only via the 
class function L x in (8) we may use the left-right in- 
variant Haar measure in (9) instead of the reduced 
Haar measure. The enhanced symmetry of the measure 
yields the following geometrical Schwinger-Dyson equa- 
tions [15], 

off. conf's {X} = J $>L cxp(-S PL ) [3L X G - (1 - L 2 X )(G' X - GS' Ph J] 



N a 



(6) 



o=l 



with Z(2) symmetric operators S a and coupling param- 
eters A a to be determined from the Schwinger-Dyson 
equations. To derive these, we proceed as follows. 

On a lattice with spacing a and temporal extent Nt 
(hence temperature T = 1/N t a) the Polyakov line is 



(11) 

Here, G[L] represents some set of functions of the 
Polyakov loop to be chosen appropriately (see below) . In 
addition, we have defined the derivative G x = dG/dL x 
and analogously for Spl- 

Switching to expectation value notation, (11) can be 
rewritten as a linear system for the couplings in (6), 

]T((1 - Ll)GS' a JX a = ((1 - Ll)G' x ) 3(L X G) . (12) 



The coefficients of this system are expectation values 
which are calculated in the full Yang-Mills ensemble ob- 
tained by MC simulation based on the SU (2) Wilson ac- 
tion. Numerically, it is of advantage to have more equa- 
tions than unknown couplings A . This is achieved by 
choosing G out of the following set of local functions, 



Ga,y € {OSa/Ly 



1. 



,N a } 



(13) 
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which represent the operators present in the equation 
of motion for L y . For fixed y, (12) then yields as many 
equations as there are couplings, namely N a . These equa- 
tions relate different two-point functions labeled by lat- 
tice sites x and y of distance r. Additional relations are 
obtained by letting the distance r run through (half of) 
the spatial lattice extent, r — 1, . . . , N s /2. Altogether, 
the overdetermined system (12) consists of N a x N s /2 
equations which are solved by least-square methods. For 
more details the reader is referred to Appendix B and 
[15]. 



III. 



CHARACTER EXPANSION 



To find a reasonable choice of operators for the ansatz 
(6) we use the beautiful analytical results of Billo et al. 
[16]. These authors have evaluated the integral (5) for the 
case at hand [S[U] being the Wilson action, X = L the 
Polyakov loop) by combining the strong coupling with 
a character expansion. For the benefit of the reader we 
briefly recapitulate their approach before we adopt it for 
our purposes. 

Recall that a character is the trace of a group element 
in an irreducible representation. If j denotes the spin of 
an SU(2) representation such that p = 2j is the length 
of the corresponding Young tableau, then the associated 
character is 



Xp (U) = tr p (U) 



sin 



((p+l)0) 



sin 6 



p = 0,1,2,... . (14) 



The characters can be entirely expressed as orthogonal 
polynomials in the traced loop L, namely the Chebyshcv 
polynomials of the second kind [17], 

This representation manifestly shows that Xp IS a P or y~ 
nomial in L of order p. The first few characters are 



Xo 



1 



X2 =4L 2 -1 



X3 



Xi = 2L , 
8L 3 - 4L 



(15) 



To determine the PLA Spl (6) in the strong coupling 
limit of the underlying SU(2) gauge theory Billo et al. 
[16] allowed for different couplings in temporal and spa- 
tial directions, denoted fit and (3 S , respectively. In terms 
of these couplings the original Wilson coupling becomes 



The formula for the temperature, 




shows that the high-temperature limit (for (3 t fixed) cor- 
responds to N t or /3 a being small. An expansion in terms 
of f3 s results in the PLA [16] 

S PL ee J2 ^ 1 1 + £ ^ P X P (L x ) Xp (L y ) 1 + . . . (16) 



P =i 



where all orders in f3 t have been summed up in terms of 
coupling coefficients n p — K p (f3 t ). The terms not writ- 
ten explicitly contain higher orders in (3 S and interac- 
tions of characters of plaquette type. The leading-order 
action (16) involves only nearest-neighbor (NN) interac- 
tions with the couplings given explicitly by 



(17) 



Asymptotically, for small (3 t , this is 
Kp = -c p p? Nt +0(p? Nt+2 ), Cp EE[2f(p+l)!]-^. (18) 

This concludes our brief discussion of [16]. In order to 
make the operator (i.e. character) content of the action 
(16) more explicit we expand the In in (16) in powers 
of [3t- From the small-/?* behavior (18), we infer that a 
product of n k's behaves as 



«P1 ■■■ K Pn =0(/3f JV *) , P =^2 Pl . 



(19) 



Thus we reshuffle the expansion of (16) such that, for 
fixed p, we first sum over all partitions of the integer p, 
then increase p by one unit, sum again etc. up to some 
maximal value, say p = 3. In this way we obtain 



Spl = S (1) + S (2) +SW + 0(f3t Nt ), 



(20) 



where is 0(f3f Nt ). Accordingly, we have a hierarchy 
of actions that become more and more suppressed 
(for small (3 t ) as p increases. We thus refer to as 
being of leading order (LO), of next -to-leading or- 
der (NLO) and so on. Abbreviating Xp x — Xp(Ae) the 
actions read explicitly 

S(1) = ^XlxXly i 

{xy) 

5(2) = ^ {n 2 XixX2y - 2 K ixLxi y ) , 

5 (3) 



(xy) 



( K 3X3xX3y - K 1 K 2 Xl X X2 x XlyX2y 
(xy) 

+ l4xLxly) ■ 

The product of characters at the same site may be further 
reduced by the 5/7(2) 'reduction formula', 

XpiX P 2 = Xpi+P2 ' Xpi+P2-2 + • • • + X\pi— P2| • 
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Note that our conventions are such that the subscripts 
get reduced by two units from left to right. Using this 
formula in the above expressions we end up with 



^ An XixXiy > 
(xy) 



(xy) 

E 

{xy) 



^ 22 X2xX2y + 2d^2X { ^X2 X , 



\(3) , \(3) 

An XlxXlw + A33 X3*X3y 



+ A13 (XlxX3y + X3xXly) 



(21) 
(22) 



(23) 



The new couplings X pq are combinations of the k's, 
namely 



A (1) 

An - 

x(2) 

A 22 - 

X (3) ... 

A 33 - 

X (3) ... 

A 13 — 



«2 
^3 



1 > 



A (2) - V 

A 02 — > 



-Kl«2 + • 



(24) 
(25) 
(26) 
(27) 



In the asymptotic regime, ft « 1, the A's may be ex- 
panded with the help of (18) (see Section V below). The 
results (21-23) from combining character and strong cou- 
pling expansion suggest the following ansatz for the PLA, 



"'i 'L T - : " E \>qXpxXqy = y ]Spq ■ ( 28 ) 
{xy),{pq) (pq) 



^pqXpxXqy — / . S p 



The couplings X pq are symmetric with respect to their 
indices pq, and p — q is even. The ansatz (28) coincides 
with the one suggested by Dumitru et al. [14] which was 
entirely based on center symmetry. It is obvious that the 
action (28) is center-symmetric as the \p are even/odd 
functions for p even/odd, x p {—L) = (—l) p x P {L). 

The terms which product pq = are localized at single 
sites and correspond to 'potentials'. Hence the action 
splits into hopping terms T and potential terms V in 
accordance with (3), 



5PL 



(pq) 

pq^O 



2YS P0 = T + V 



p 



(29) 



With Xo = 1 the potential V has the explicit form 

V = 2d Aop^Txp* • (30) 

p even x 

For what follows we need some notation. We will, of 
course, truncate our actions at some maximum 'spin' p, 
say at p = r. Thus we define the truncated actions 

r r r 

S r = S pq = S pq + 2 S p o = T r + V r , (31) 

(pq) (pq) p=2 

pq^O 



where all summations are cut off at p = r. Explicitly, the 
first few terms are 

51 = Sn , (32) 

5 2 = (Su + S 22 ) + 2S 20 = T 2 + V 2 , (33) 

5 3 - (T 2 + S31 + S 33 ) + 2S 20 = T 3 + V 3 . (34) 

Note that actually V 3 = V 2 . In the strong-coupling ex- 
pansion the action S r is of order 0([3 rNt ). As before, we 
refer to p = 1 (fundamental representation) as the LO, 
to p = 2 (adjoint representation) as NLO and so on. In 
this paper we will not go beyond p — 3, a truncation that 
neglects terms which are NNNLO in the strong-coupling 
expansion. 

According to (33) the first potential term arising in the 
strong-coupling character expansion is of NLO (f3 2Nt ) 
and quadratic in the Polyakov loop L, 



V 2 = 2d\ a2 Y X2a 



(35) 



IV. MEAN FIELD APPROXIMATION 

Before we actually relate the PLAs (28) to Yang-Mills 
theory let us analyze their critical behavior which is in- 
teresting in itself. Both via mean-field (MF) analysis 
and MC simulation we will see that the models typically 
have a second order phase transition at certain critical 
couplings. Obviously this should match with the Yang 
Mills critical behavior. 

The models involving more than a single hopping term 
also show a first order transition at which the order pa- 
rameter (L) jumps. For the SU{2) case discussed here 
this transition is not related to Yang-Mills. Matching 
to the latter hence implies that the effective couplings 
should stay away from the first-order critical surface. 

To develop a MF approximation for the Polyakov loop 
L with its nontrivial target space [— 1 , 1] we use a varia- 
tional approach based on the text [18]. Our starting point 
is the effective action Y for the Polyakov loop dynamics. 
This time, the term 'effective action' refers to the gener- 
ating functional for the 1PI correlators of the Polyakov 
loop. The former represents the complete information of 
the quantum field theories based on the PLAs Spl- This 
is obvious from the fact that T is the Legendre transform 
of the Schwinger functional W, 

T[L] = (&W)[L] , W[J]=logZ[J], (36) 
Z[J] - J @L exp ( - S PL [L] + (J, L)) . (37) 

These are the standard definitions to be found in any text 
book on quantum field theory. For an alternative, varia- 
tional characterization of F [18] we consider the following 
probability measures on field space, 



dn[L] = ®Lp[L] , 



(38) 
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with nonnegative function p[L] and S>L as in (10). Aver- 
ages are calculated with //. For example, the mean PLA 
is 

(S PL )» = J dv[L}Sp L [L} , (39) 

while the Boltzmann-Gibbs-Shannon entropy is given by 
the average of logp, 

SbgsM = -(logP>M = " / ML] \ogp[L} . (40) 

The relevant variational principles arc obtained as fol- 
lows. By subtracting (40) from (39) one forms 

^V] = (Sp L +logp) M = F[p]. (41) 

This analog of the free energy is varied with respect to p 
under appropriately chosen constraints. These are added 
via Lagrange multipliers. If one just requires normaliza- 
tion to unity, (1) M = 1, one finds that the probability p 
for which (41) becomes extremal is given by the standard 
measure, p = exp(— Spl)/Z[0]. Inserting this into (41) 
yields the infimum 

mfF\p] = -logZ[0] = -W[0] . (42) 

Comparing with (36) this may be interpreted as the effec- 
tive action for J = L = 0. If we vary F in (41) keeping 
the expectation value of L fixed at L x by means of a 
Lagrange multiplier J x , we find the probability 

p -S+(J,L) 

P[L] = -^p , (43) 

where J is to be viewed function of L, obtained 

via inverting the implicit relation L = 5W/SJ. Plugging 
(43) into (41) yields a new variational infimum which is 
precisely the effective action, 

inf {F[fj] | {L x )n = L X } = (J, L) - W[J] = T[L] . (44) 

In the MF approximation for the effective action one min- 
imizes only with respect to all product measures v€p, 

dv\L\ = Y\_dv x (L x ) , dv x {u) = p x (u) dh(u) , (45) 

X 

where h is the reduced Haar measure from (10). Hence, 
the exact expression (44) is replaced by its MF version 
according to 

mf (F[u] | {L x ) v = L x ) = Umf[L] . (46) 

Clearly, from the variational principle, the MF effective 
action Umf bounds the effective action T from above. 
Since the set of all product measures is not convex (unlike 
the set of all probability measures), the MF action need 
not be convex. We may, however, use its convex hull 



given by the double Legendre transformation, Fmf = 
-5? 2 (Umf), which for non-convex Umf will be a better 
approximation, T < Fmf < ^mf- 

For a product measure the entropy and mean action 
entering F in (46) turn into sums (of products) of single 
site expectation values. With the abbreviation 

J du x {u)f{u) = (f) x , (47) 

we find for the general class of character actions (28) 

(xy),{pq) x 

This becomes extremal for the single site measure 

*.(«) = ^IZX^+M , (49) 

Z{){J X ) 

which replaces (43). V denotes the potential from (30), 
and the normalization factor is the single-site partition 
function, 

zo(j) = J dh(u) cxp (-V(u)+ju) = cxp (w (j)) . (50) 

In analogy with (43) the Lagrange multiplier (or exter- 
nal source) J x in the single-site measure p x has to be 
eliminated. This is done by inverting the relation 

L x = w' (J x ), (51) 

(the prime henceforth denoting derivatives with respect 
to J) so that Jx = J X (L X ). Since the Schwinger function 
is strictly convex, the relation between the mean field 
L and the external source J in (51) is one-to-one. This 
will become important in a moment. 

Taking these considerations into account the infimum 
of (48) becomes the MF potential 

U M f[L} = X pq(x P )x(Xq)y + ^2lo(L x ), (52) 

(xy) x 

(pq)^O 

where all single-site expectation values are calculated 
with the measure (49) subject to the condition (51). The 
quantity 70 is the Legendre-transform of Wo, 

70 W = mf (tj - w (j)) ■ (53) 

In this paper we consider effective potentials rather than 
effective actions and hence assume that the source and 
mean field are constant, L x = I. The effective potential 
is just the effective action for constant fields, divided by 
the number of lattice points. The MF potential (52) then 
simplifies to 

umfW = d \ q (x P )(x q )+lo(e) ■ (54) 

(P9)#0 



6 



This bounds the true effective potential from above and 
so does its convex hull, 



7mf = -S? 2 (wmf) . 



(55) 



In [19] it is proved that umf is the MF approximation to 
the constraint effective potential [20]. To calculate umf 
one proceeds as follows: 

1. For a given source j one computes the single-site 
partition function zo according to (50) and its log- 
arithm, the Schwinger function, wq = logzo- 

2. The minimizing j in (53) satisfies the equation 

which must be inverted to yield j — j{tj. 

3. The solution is used to calculate 70, 

10 (e)=j(e)£- Wo (j(e)). (57) 

4. Then one computes the expectation values 

(Xp) = J exp (ju-V(u))x P (u) dh(u) , (58) 

where j = j(£) as obtained from (56) is inserted 
everywhere. Plugging (58) into (54) finally results 
in the MF potential mmf- 

We are interested in the MF expectation value of the 
Polyakov loop which minimizes the MF potential (54). 
Since the relation I = £(j) in (56) is one-to-one, the 
condition that £ is a minimum of Umf is equivalent to 
u 'mf = 0- We thus need the j-derivatives of 70, 



1 (new d lod£ j , 
^ i{ ^ = ^£dj=2 {Xl) ' 



(59) 



as well as of the mean characters, 

2<XpY = (Xp+i) + <X P -i> - (XpXXi) • (60) 
Setting u' MF = results in the self-consistency condition 

J2 2rfA M (( Xp )(x g )) / +i(Xi)' = 0, (61) 

(pq)^O 

where all j-derivatives are calculated via (60). 

Let us finally find the couplings for which the curvature 
of umf at the origin changes sign. For these, the value 
of «mf at £ = turns from a minimum to a local maxi- 
mum signaling a second order phase transition. Since the 
potential in (49) is center symmetric, the strictly convex 
functions w and 70 are both even functions with ab- 
solute minimum at the origin. It follows that £ = is 
mapped to j = and that (for these arguments) the 



second ^-derivative of mmf is proportional to its second 
j-derivative, 



d 2 u M F 



d£ 2 







d 2 U M F 



3=0 



= 



(62) 



To calculate the second j-derivative of umf we need 

4 (x P >"Uo = ((Xp+2+X P +X P -2) - (Xp)<X2))| j=0 , (63) 

which holds for even potentials. Together with (60) and 
(61) this implies the 'zero-curvature condition' 

J2 2dA pg ((x P )(x 9 })"b=o + (xi)1 j=0 = 0. (64) 

(p 9 )#0 

A. Ising-type models 

If there is no potential in the PLA (29) and if the 
hopping terms contain only NN interactions we refer to 
the resulting actions as Ising-type models. The associated 
actions have been denoted T r in (31). Setting V = in 
(50) we obtain 

zo(j) = Io(j) - h{j) , (65) 
while the mean characters in (58) become 



(y ) = A (7) - U?)- WJ) 



(66) 



For p —1 we find the field conjugate to the source j, 

h{j)-h{j) 



(Xi) = 2£(j) = 



(67) 



IoU)-h(j) ' 

The inverse relation j = j(£) is used in the MF potential, 
umf(£) =dJ2 \ q A p(j)\(J) + V - Mj) , (68) 

pq^O 

the minimum of which is determined by the self- 
consistency condition (61). Since A p ^ (0) and Ap_^ 1 (0) 
both vanish, the condition (64) for a second order phase 
transition simplifies to 



(69) 



Hence, in the MF approximation all Ising-type models 
show a second order phase transition at the critical cou- 
pling Aii c . In the presence of NLO couplings the possi- 
bility of first-order phase transitions arises. In this case, 
the critical couplings have to be determined numerically 
(see below). 

Historically, the first derivation of an effective action 
for the confinement -deconfincment phase transition is 
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due to Polonyi and Szlachanyi [21]. They have also uti- 
lized the strong-coupling expansion on a Euclidean lat- 
tice and already found the LO contribution from (32), 



Si = T\ = y Anxi x xiy 



(70) 



Note that our sign-convention is such that negative An 
corresponds to 'ferromagnetism'. The action (70) entails 
the MF potential 



u ( ^ F (e) = 4d\ n e 2 + e 3 - w a (j) , 



(71) 



where one inverts the map j i— ► I in (67) to obtain j(£). 
This potential is minimal for £ — L which solves 



8d\ n L + j(L) = , 
implying the self-consistency condition 
2i = Ai(-8dAuZ) . 



(72) 



(73) 



The effective potential and order parameter are plotted 
in FIG. 2. Near the critical coupling the order parameter 




FIG. 2: MF prediction for the effective potential and order 
parameter for the Ising-type PLA Si = T\ from (70). 



has the typical square root behavior, 
L ? 



An 



An 



1 for |An|\|Aiic|. 



(74) 



Apart from the second-order transition at An c = — l/2d 
there is no further phase transition (in the MF approx- 
imation) for the LO ansatz (70) involving only one cou- 
pling. 



As already mentioned, the situation is different for the 
Ising-type PLA with NLO coupling A 2 2, 



T 2=2_j (*>ilXlxXly + A 2 2X2xX2y) 

The corresponding MF potential reads 



r,( 2 ) 

U MF 



AdXni 2 + dX 22 A 2 2 (j) + £j - woU) 



(75) 



(76) 



where A 2 is defined in (66) and j(£) follows from (67). In 
addition to the second-order phase transition at An c = 
— 1/6 the model has a first order transition: the order pa- 
rameter jumps as a function of A 22 as long as An > An c . 
The behavior of the order parameter as a function of the 
two couplings is shown in FIG. 3. 




FIG. 3: MF behavior of the order parameter L for the Ising- 
type PLA T 2 from (75). 

In FIG. 4 the results for the order parameter L as ob- 
tained by MC simulations are plotted. Shown is L for sev- 
eral thousand sample points in the (An, A 22 ) plane. Near 
the second-order phase transition curve with An < 0.2 
about 2.5M updates were sufficient, whereas near the 
transition curve with An > 0.2 at least 50M updates 
were required. 




FIG. 4: MC data for the order parameter L for the Ising-type 
PLA Ti from (75). Shown are the critical curves of the MF 
approximation (dark) and MC simulations (bright). 

We interpret this slow convergence as a signal for a first- 
order phase transition as predicted by the MF approxi- 
mation, see FIG. 4. 
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B. Ginzburg— Landau— type models 

Inspired by the nomenclature in [9] we call effective ac- 
tions with a nonvanishing potential as in (29) Ginzburg- 
Landau-type models. The lowest-order potential term is 
V 2 given in (35) which is of NLO in the strong-coupling 
expansion. Hence we first consider the action 



S 2 =T 2 + V 2 , V 2 = 2d\ 02 X2x • 



(77) 



For A 2 7^ the MF potentials and mean values must 
be computed numerically. Only for the special case of 
vanishing source may the relevant integrals be found an- 
alytically, 



dh(u) e- 2au ~ X 2p(u) = [J p (a) + I P+ i(a)} . (78) 



For what follows it is useful to define 

Spq = Ip(4d\o 2 ) + Iq{4:dXo 2 ) , 

d pq = I p (4d\ Q2 ) - Iq(4d\ 02 ) , 
so that the mean characters can be written as 

<x 2p >U = (-F^f. 



(79) 
(80) 



(81) 



Using (79) and (80) the critical surface for second-order 
transitions in the 3-dimensional space of coupling con- 
stants An, A22, A 2 is determined by the equation 



An =- 



1 Soi 



+ A 2 



2d d 02 

which is plotted in FIG. 5 



(£) 2 


f £03 


S\ 2 \ 




\S\ 2 


soi) 



(82) 




FIG. 5: Critical surface for the Ginzburg-Landau PLA S2 
from (77). 

The NNLO Ginzburg-Landau-type model is given by 
S3 = T a + V 3 , (83) 



and has the same potential V3 = V 2 as the NLO action 
(77) whence (81) still holds. It follows that the critical 
surface for second-order transitions in the 5-dimensional 
space of coupling constants An, A22, A 2, A13, A33 is deter- 
mined by the transcendental equation 



A 
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1 got 
2d d Q2 

— A33 ( —, — 

d 02 



\ 2 2 



{d 02 ) 


( S03 


S\ 2 \ 




\S\ 2 


soi) 



2Ai 3 ^ ■ 
d 02 



(84) 



For fixed small couplings A33 and A13 this surface is very 
close to the one in FIG. 5. 

Summarizing we have seen that the effective models 
undergo both first- and second-order phase transitions, 
and that an MF prediction for the phase structure of 
these models agrees well with the results of MC simula- 
tions. 



IMC RESULTS FOR \ pq (/3) 



Employing our IMC method based on the geometrical 
Schwinger-Dyson equations of Section II we have calcu- 
lated the couplings \ pq for different PLAs as functions of 



j3. For our lattice with N s 
cal Wilson coupling is (i c = 
hierarchy 



= 20 and N t = 
2.30. Below f3 c 



I An I > |A 22 | , |A 02 | > |A 33 | , |Ai 3 | 



4 the criti- 
wc find the 



(85) 



in agreement with the strong-coupling expansion (21-23). 
According to Ogilvie [22] , the weak-coupling asymptotics 
of An is linear in /3, 



An(/3) = -- 



2N c N t 

which in our case (Nq = 2) leads to 



+ const (P > p c ) , 



An(/3) = - 



L 
16 



const 



-0.0625 3 + const 



(86) 



We have compared our IMC results for the couplings 
Xpq(P) with the strong-coupling predictions (17, 24-27) 
and the weak-coupling result (86). As expected from our 
reasoning above, the lowest order PLAs based on group 
characters approximate the true Polyakov-loop dynamics 
very well in the strong-coupling regime. For weak cou- 
pling we find the linear relation (86) already for the LO 
PLA. For the NLO Ginzburg-Landau-type model with 
3 couplings the slope —0.0614 is very close to the weak 
coupling result —0.0625 in (86). Thus we are confident 
that our PLAs describe the true Polyakov loop dynamics 
below and above the critical Wilson coupling very well. 



A. Leading— order action 

The effective couplings An for the LO Ising-type 
model (70) for /^-values below and above the critical 
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TABLE I: Comparison of the critical coupling values for the 
LO Ising-type model Si. The 'exact' value is obtained via 
MC simulation. 



0.0 



method 



critical coupling 



MC simulation 
MF 

strong coupling 
IMC 



Auc = -0.18 
Aiic = -0.17 
Auc = -0.11 
An (/3 C ) = -0.14 



-0.1 



f3 c = 2.30 are listed in TABLE VI, Appendix A. We read 
off the value An(/3 C ) = —0.132. If Si would be the ex- 
act PLA then its critical coupling An c would be An(/3 C ). 
A direct MC simulation of the action Si reveals that 
this model has critical coupling An c = —0.18. The MF 
prediction An c = —0.17 comes surprisingly close to the 
former 'true' value. The critical coupling may alterna- 
tively be estimated by using the strong-coupling results 
(18) and (24) to calculate An(/3 C ) (extrapolating them to 
N t = 4 and [3 C = 2.30), 



-0.2 - 




An = -G9/4)"' « -0.11 . (87) 

The output of the different methods is compiled in TA- 
BLE I. The values obtained are quite close to each other. 
The 'true' value —0.18 stems from simulating (70) with 
the MF approximation coming closest. Somewhat sur- 
prisingly, the strong-coupling result yields the right order 
of magnitude. The discrepancy between direct simulation 
and IMC just means that the action (70) represents an 
oversimplification and does not match with Yang-Mills 
well enough. The IMC value An(/3 C ) = —0.14 constitutes 
a compromise equivalent to a one-parameter fit of the ef- 
fective to the Yang-Mills Schwinger-Dyson equations or 
n-point functions. 

In FIG. 6 we compare the results for An from strong- 
and weak-coupling expansion with the data of our IMC 
simulations. Note that the asymptotic formula (87) 
works up to f3 ~ 2. The IMC values for the LO action 
Si in (70) are marked with circles. Note further that the 
dependence An(/3) is indeed linear in the weak-coupling 
regime, An = —0.0392 f3 — 0.0774, in accordance with the 
prediction (86) for the weak-coupling asymptotics. The 
slope, however, turns out being too small. This will be 
remedied in what follows by including more couplings. 

B. Next— to— leading— order actions 

The effective action (70) has been confirmed and ex- 
tended by several authors [22-25]. We have also checked 
its generalizations by considering the NLO Ising-type 
model (75) without potential terms and the Ginzburg- 
Landau action (77) containing all NLO terms from (33). 

The NLO IMC results for An(/3) are displayed in 
FIG. 6 (crosses and asterisks), those for \22(/3) in FIG. 7 



FIG. 6: The coupling An as a function of (3. Dashed curve 
for (3 < p c : asymptotic behavior (87); dashed lines for /3 > f3 c : 
linear fits A = af3+b; data points: IMC results from Tables VI- 
VIII in Appendix A. 



TABLE II: Slopes of the linear fits to the weak-coupling 
asymptotics (86) of An. 



model 


l/2N N t 


5*i 


T 2 


s 2 


slope 


-0.0625 


-0.0392 


-0.0505 


-0.0614 



(same symbols) and those for Xq2(P) in FIG. 8 (aster- 
isks). The error bars for the couplings are listed in the 
tables in Appendix A. Upon comparing the values for An 
in FIG. 6 (or TABLE VI) we note that for < (3 C the 
predictions for An are almost model independent. In the 
weak-coupling regime, on the other hand, they are less 
stable. Hence, adding terms to the PLA may change the 
coupling constants considerably in this regime. 

For the NLO actions the dependence Xn(f3) is linear in 
the weak coupling regime, similarly as for the LO action 
Si. However, the slopes a in the linear fits An (j3) = af3+b 
above (3 C are model dependent. For the LO and NLO 
PLA they are given in TABLE II together with the 
weak-coupling result (86). For the Ising-type models 
without potential terms the slope is not reproduced very 
well, and we conclude that in the weak-coupling or high- 
temperature regime potential terms should be included 
for an accurate description of the Polyakov-loop dynam- 
ics. Indeed, the slope for the action (77) with potential 
term V2 is almost identical to the prediction (86) of the 
weak-coupling asymptotics. 

FIG. 7 shows the dependence of the adjoint coupling 



constant A22 on the Wilson coupling j3. Also shown is the 
prediction (25) of the strong coupling expansion using 
(18), 

A22(/3) --ii(i) 8 "-°- 69753 (0 8 ' (88) 

which again reproduces the IMC results for /3 < 2. 
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FIG. 7: The adjoint coupling A22 as a function of /3 for the 
JVLO and NNLO actions. Dashed curve: asymptotic behavior 
(88); data points: IMC results from Tables VI-VIII. 

The critical coupling 

A 22 (/3 C ) = -0.008[1] (89) 

is the same for both PLAs (70) and (77). Up to (3 C the 
two actions have almost identical couplings A22 ■ It seems 
that for (3 > (3 C these couplings also depend linearly on 
/?, as was the case for An(/3). 

C. Next-to-next-to- leading order 

We have seen that the NLO approximations describe 
the Polyakov loop dynamics very well in the sym- 
metric strong-coupling phase. In the weak-coupling 
regime, however, there is still room for improvement. 
Hence we have calculated the five coupling constants 
An, A22, A02, A13, A33 appearing in the general NNLO 
PLA S3 for several values of the Wilson coupling. This 
action is the sum of all terms up to order (3f Nt in the 
strong-coupling expansion. As expected, adding the 
third order terms S31 and S33 does not change the lower- 
order couplings An, ^22j A02 (as obtained via S2) in the 
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FIG. 8: The coupling A02 as a function of the Wilson cou- 
pling (5 for NLO, NNLO and Ginzburg-Landau actions; data 
points: IMC results from Tables VII-IX. 

broken phase. This can be seen in FIG.s 6, 7 and 8, 
where the IMC results for the NNLO action S3 from (34) 
are depicted by triangles. The numerical values for these 
couplings and the couplings A33 and A13 together with 
their statistical errors are given in Appendix A. 

D. Ginzburg-Landau models 

In his review [9] , Svetitsky has suggested to emphasize 
the potential term in (3) by specializing to an ansatz of 
Ginzburg-Landau type, 

Si + . ( 9 °) 

X 

with center-symmetric potential (30). Replacing the co- 
efficients Ao p of the characters in this potential by 

A 4 = 2d ■ 4 2 A 04 , (91) 
A 2 = 2d ■ 4(A 02 - 3A 04 ) , (92) 

leads for p < 4 to the even polynomial 

V(L) = \ 2 L 2 + A4L 4 + const . (93) 

In Table IX we have listed the couplings for the models 
with quadratic and quartic potentials, 

Sg) = S 1 + X 2 J2Ll , (94) 

X 

x 
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TABLE III: Effective couplings of different PLAs at /3 = /3 C = 
2.30. 



model 


Si 


T 2 


s 2 


s 3 


a(2) 
GL 


An 03 c ) 


-0.133 


-0.135 


-0.156 


-0.156 


— 0.155 


A 22 (/3 C ) 




0.006 


0.006 


0.006 




A o2 03 c ) 






0.011 


0.012 


0.012 



obtained via IMC within our standard range of j3. The 
values for An both for strong and weak coupling are al- 
most identical to those of the NLO model 82- For this 
reason we have refrained from plotting them in FIG. 6. 
The potential couplings are important in the broken 
weak-coupling phase where they become sizable. The 
coupling A02 for the Ginzburg-Landau models (94,95) is 
shown in FIG. 8. 



E. Summary 

The values for the couplings of the different PLAs aris- 
ing at critical Wilson coupling (3 C = 2.30 are listed in 
TABLE III. They are almost model independent. An, 
in particular, is always close to the MF value —0.167 (for 
the Ising type models). 

The couplings below and above (3 C and their statisti- 
cal errors are compiled in Appendix A. There one may 
also find A33,Ai3 and A04 = A 4 /32d for different Wilson 
couplings. 

The stability of the couplings for (3 < (3 C is a strong 
indication that (in this regime) the Yang-Mills ensemble 
is very well approximated already by the NLO models 
with 2 or 3 couplings. The results of the following section 
will further confirm this statement. 

The Ising-type coupling An becomes a linear function 
of (3 in the weak-coupling regime, in accordance with 
the weak-coupling prediction (86). For the NLO ac- 
tion the slope is —0.0614 which compares favorably with 
the weak-coupling value —0.0625. For the Ginzburg- 
Landau-type actions the slope is almost identical to the 
one of the NLO models. 

The Ising-type couplings change rapidly at the critical 
Wilson coupling j3 c — 2.30 as demonstrated in FIG.s 6 
and 7. For example, the coupling An decreases from —0.1 
below (3 C to —0.2 above (3 C . This jump of An forces the 
system into the ferromagnetic phase. For A22 the jump 
is even more dramatic, from to 0.04. The potential 
couplings A02 and A04 change more smoothly when the 
systems changes from the symmetric to the broken phase. 



VI. TWO-POINT FUNCTIONS 

Let us finally check the quality of our PLAs which, 
after all, should represent approximations to Yang-Mills 



TABLE IV: Effective couplin; 


5s for NLO actions at (3 = 2.2. 


P = 2.2 




An 


A22 A20 






0.1033 


0.0019 


s 2 




0.1168 


0.00042 0.0064 


TABLE V: Effective couplings 


for NLO and NNLO actions at 


P = 3.0. 








p = 3 


An 


A22 


A 2 o A33 A31 


T 2 


0.2207 


0.0320 




s 2 


0.2361 


0.0194 


0.0165 


S3 


0.2506 


0.0311 


0.0236 0.0035 0.0040 



theory. To this end we compare the Yang-Mills two- 
point function at different /9-values with those of our 
effective models inserting the couplings X pq (f3) obtained 
via IMC. For the NLO actions at j3 = 2.2 these couplings 
are displayed in TABLE IV, while for (3 = 3.0 we have 
the NLO and NNLO couplings of TABLE V. With these 
couplings we have simulated the models with actions 



T 2 



So 



'GL i 



S, 



(2) 
GL 1 



(96) 



and calculated the two-point functions displayed in 
FIG.s 9 and 10. As expected, the agreement in the 
center-symmetric phase ((3 — 2.2) is very good, while 
deep in the broken phase (/3 = 3.0) there appears to be 
room for improvement. For (3 — 2.2 the two-point func- 
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FIG. 9: The Yang-Mills (YM) two-point function compared 
to the ones obtained from the NLO effective actions T 2 and 
S 2 ; (N B ,N t ) = (20,4). 



tions of the three effective models considered are almost 
identical to the Yang-Mills two-point function. The data 

(2) 

points for S3 and Sq^ cannot be distinguished from those 
for the NLO model S2 and hence are not displayed in 
FIG. 9. 

For (3 = 3.0 the two-point functions and the expec- 
tation value of the mean field are model dependent. In 
FIG. 10 we have plotted the two-point function of Yang- 
Mills theory and of the NLO and NNLO effective actions. 

For the NLO approximation T 2 the value for the con- 
densate is approximately 20 percent below the Yang- 
Mills value. Including potential terms (5*2) and NNLO 
terms (S3) improves the approximation somewhat as 
FIG. 10 shows. The two-point function for the 
Ginzburg-Landau action S G -[ is almost identical to the 
one of T2 (and thus has not been displayed). This im- 
plies that also higher-order Ising (or hopping) terms are 
important when (3 > (3 C . 
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FIG. 10: The Yang-Mills (YM) two-point function compared 
to the one obtained from the NLO (T2 and S2) and NNLO 
(S3) effective actions; (N s ,N t ) = (20,4). 



VII. DISCUSSION 
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characters, the couplings have been determined by our 
standard inverse Monte-Carlo routines. We found that 
the couplings decrease rapidly not only if we go to higher 
representations (i.e. larger p) as above, but also if we in- 
crease the number of links within the plaquette operators 
used. The leading r = 1 term Si of (32) with coupling 
An dominates by one order of magnitude compared to 
the terms with r — v2- This clearly indicates that the 
effective interactions are short-ranged in accordance with 
the Svetitsky-Yaffe conjecture. Simulating the effective 
action including plaquette operators with the couplings 
obtained via IMC we have calculated the two-point func- 
tion in both phases. Including a total of 14 couplings the 
matching between Yang-Mills and the effective action be- 
comes perfect both in the broken and symmetric phase 
[15]. It should, however, be stressed that this brute- 
force numerical calculation does not provide too much of 
physical intuition. We believe that the present paper, in 
particular the mean-field analysis, improves upon [15] in 
this respect. 

The most straightforward generalization of our analy- 
sis obviously is to go to SU (3) Yang-Mills theory where 
one expects a first order phase transition. Work in this 
direction is under way. 



APPENDIX A: TABLES 

The following TABLEs VI-IX contain the effective 
couplings X pq with statistical errors for various values of 
(3. 



TABLE VI: IMC results for the effective coupling An in the 
LO action Si and the two couplings An, A22 in T 2 . 



Yang-Mills- f3 


Si : An 


Ti : An 


T2 : A22 


1.70 


-0.0305[4] 


-0.030[2] 


-0.001[1] 


2.20 


-0.1040[2] 


-0.103[1] 


-0.0019[7] 


2.28 


-0.1231[3] 


-0.123[1] 


0.000[1] 


2.29 


-0.1267[4] 


-0.127[1] 


0.002[1] 


2.30 


-0.1325[5] 


0.135 [1] 


0.006[1] 


2.32 


-0.1512[3] 


-0.1598[9] 


0.023[1] 


2.34 


-0.1585[2] 


-0.1697[6] 


0.030[1] 


2.38 


-0.1658[1] 


-0.1799[4] 


0.037[1] 


2.60 


-0.1792[1] 


-0.1989[4] 


0.0382[8] 


2.80 


-0.1874[1] 


-0.2103[5] 


0.0350[8] 


3.00 


-0.1952[1] 


-0.2206[6] 


0.0319[8] 


3.50 


-0.2146[1] 


-0.2446[6] 


0.0264[6] 


4.00 


-0.2343[1] 


-0.2670[9] 


0.0225[6] 



To really obtain satisfactory approximations to Yang- 
Mills expectation values (e.g. two-point functions) for 
all /3-values one has to go beyond nearest-neighbor in- 
teractions in the effective theory. This has been done in 
[15] where operators with r = \x — y\ up to V2 ('pla- 
quette operators') were included. We briefly recapitulate 
the results of this brute-force calculation. By reformu- 
lating the Schwinger-Dyson equations (12) in terms of 



APPENDIX B: NUMERICAL DETAILS 

All MC calculations have been performed on a 20 3 x 
4 lattice for which the critical Wilson coupling is (3 C — 
2.30. The simulations have been done for (3 ranging from 
1.1 to 4.0. We have used a standard 'pseudo-heatbath' 
algorithm [26, 27] due to Miller [28]. 
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TABLE VII: IMC results for the three couplings in S 2 . 



p 




An 




A22 


A20 


1.10 


-0 


004 [4] 


-0.000(2] 


-0.0012(5] 


1.70 


-0 


030 [4] 


-0.001 2 

L J 


-0.0007(5] 


2.20 


-0 


116(2] 


0.000(2] 


0.0064(5] 


2.28 


-0 


144(2] 


0.002(2] 


0.0101(5] 


2.29 


-0 


148 [2] 


0.004(2] 


0.0108(5] 


2.30 


-0 


156(2] 


0.006(2] 


0.0114(5] 


2.32 


-0 


178(2] 


0.017(2] 


0.0125(6] 


2.34 


-0 


187(1] 


0.022(2] 


0.0130(6] 


2.38 


-0 


195(1] 


0.026(2] 


0.0130(6] 


2.60 


-0 


212(1] 


0.026(2] 


0.013(1] 


2.80 


-0 


224(1] 


0.022(2] 


0.014(1] 


3.00 


-0 


236 [1] 


0.019(2] 


0.016(1] 


3.50 


-0 


267(1] 


0.010(2] 


0.025(1] 


4.00 


-0 


299(1] 


0.003(2] 


0.037(1] 


TABLE VIII: IMC results for five couplings 


m 03. 


P 


An 


A22 


A20 


A33 


A31 


1.10 


-0.004[1] 


0.000[1] 


-0.0011(2] 


0.000(1] 


0.0003(6] 


1.70 


-0.029(2] 


0.000(1] 


-0.0006(2] 


0.0001(8] 


0.0000(6] 


2.20 


-0.116(1] 


0.000[1] 


0.0064(2] 


0.0015(7] 


0.0000(6] 


2.28 


-0.143(1] 


0.002(1] 


0.0102(3] 


0.0015(8] 


0.0002(6] 


2.29 


-0.148(1] 


0.003(1] 


0.0108(3] 


0.0015(8] 


0.0000(6] 


2.30 


-0.156(1] 


0.006(1] 


0.0115(3] 


0.001(1] 


0.0000(6] 


2.32 


-0.179(1] 


0.019(1] 


0.0129(3] 


-0.0006(9] 


-0.0004(5] 


2.34 


-0.1892[8] 


0.026(1] 


0.0138(3] 


-0.0014(8] 


-0.0014(6] 


2.38 


-0.1999[6] 


0.035(1] 


0.0148(4] 


-0.0032(9] 


-0.0026(5] 


2.60 


-0.2227(8] 


0.040(1] 


0.0176(5] 


-0.005(1] 


-0.0041(4] 


2.80 


-0.236(1] 


0.035(1] 


0.0202(8] 


-0.0041(9] 


-0.0040(4] 


3.00 


-0.250(1] 


0.031(1] 


0.0235(8] 


-0.0034(8] 


-0.0039(4] 


3.50 


-0.283(1] 


0.020(1] 


0.032(1] 


-0.002(1] 


-0.0032(6] 


4.00 


-0.291(3] 


0.003(1] 


0.029(2] 


-0.0033(9] 


0.0029(8] 



The IMC routine has been implemented as follows. For 
each action term S pq and site we have chosen the operator 



which leads to the Schwinger-Dyson equations 



(Bl) 



^pq(0-~ Lx)Gpq,zS p q t x) — ((1 — L x )G tX — 'iL x G) . 

(pq) 

(B2) 

Due to translational invariancc the coefficients of G z and 
G z i are equal if \x — z\ = \x — z'\. In order to have a suf- 



ficiently overdetermined system (for fixed pq) we choose 
the N s operators G pq ^, d = 0...N S . Independent of 
our choice of PLA we have always used the full set of 
operators up to truncation values p, q = 5, i.e. 



Gu , G22 , G20 , G33 , G31 , G44 , G42 , G40 , G55 



(B3) 



with < d < 8 leading to a total of 81 operators (and 
equations). Translational invariance admits to use the 
spatial average of each Schwinger-Dyson equation and 
every configuration. The overdetermined system is then 



TABLE IX: IMC results for the couplings Aii,A 2 ,A 4 in S, 



(2) 



2.20 
2.25 
2.28 
2.29 
2.30 
2.32 
2.34 
2.40 



A 4 



-0.110 
-0.119 
-0.127 
-0.128 
-0.157 
-0.173 
-0.176 
-0.190 



0.186 
0.237 
0.288 
0.303 
0.453 
0.621 
0.698 
0.979 



-0.002 
-0.003 
-0.006 
-0.007 
-0.020 
-0.057 
-0.093 
-0.256 



solved via least-square methods. We have checked that 
the couplings obtained in this way follow a normal distri- 
bution, as expected. Hence we calculated the standard 
deviation a and took la as our error. Autocorrelation 
effects have been eliminated via binning. Our statistics 
(5k to 10k configurations) entail a statistical error of 10~ 4 
which translates into an uncertainty for the couplings in 
the NLO action of the order of a few percent. The NNLO 
couplings A 33 and A31, however, have statistical errors of 
about 20%. 

Systematic errors are mainly due to the dependence 
of the couplings on the operator bases used in the 
Schwinger-Dyson equations. 
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